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Abstract. The anomalous dimensions of the general flavour non-singlet quark bilinear cur- 
rents ^^^^"^ . . . 7'^"I'i/; are computed at three loops in quantum chromodynamics in the minimal 
subtraction scheme. The dimension of the tensor current emerges for the case n = 2 and the 
anomalous dimension for the general flavour singlet current is also discussed. 
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In recent years the construction of perturbative results at large numbers of loops in quantum 
chromodynamics, (QCD), has advanced substantially. For example, the /3- function and various 
other renormalization group functions are now known to four loops in the MS scheme, [|^, ^, ^. 
One motivation for such calculations rests in essence with the improvement in precision of 
experiments relating to the strong nuclear force. Whilst such four loop results represent the 
current limit of our ability to calculate in perturbative QCD, various quantities remain to be 
determined at three loops. For instance, whilst the anomalous dimension of various quark 
bilinear currents such as the axial vector current are known to at least three loops, |5|, that 
of the closely related flavour non-singlet tensor current Tpa^'^^ where a^^ = [7^,7^^] and ip is 
the quark field, is only known at two loops. This anomalous dimension occurs in the matching 
between currents in QCD heavy quark effective theory, [^. In the two loop MS result was 
deduced indirectly as opposed to performing an explicit two loop renormalization. In this letter 
we complete the gap in the literature by providing the value of the three loop MS anomalous 
dimension of the current tpa'^'^ip. While this may appear to represent a moderate progression 
in this area, we choose to determine it not by explicit calculation of the current on its own but 
deduce it as a corollary of the renormalization of a set of generalized currents. These are j^^j"^" 
= 'tl>T^^^"^"ip, where T^^^"^" is the antisymmetric product of n 7-matrices. Clearly for QCD 
only the first five currents are non-zero but our motivation for considering the larger set relates 
to other issues. First, it is clear that for n = 2 the result we seek for the tensor current will 
emerge simply. However, by calculating with the general currents the results for n = and 1 
and the naive results for n = 3 and 4 will provide non-trivial checks on the calculation from 
the point of view of, say, symmetry factors and correctly including all the three loop Feynman 
diagrams. These three loop results have been determined in |^ ^, ^. We have refered to the 
results for n = 3 and 4 as naive since in these cases there is a connection with the 7^ problem of 
dimensional regularization. In particular the anomalous dimensions which emerge for n = 3 and 
4 do not correspond directly with the anomalous dimensions of the axial vector and pseudoscalar 
currents respectively. The correct anomalous dimensions for these currents are determined by 
including a finite renormalization to ensure that chiral symmetry is preserved, Such a 
procedure was elaborated on in detail in [|| but this issue in the context of the matrices r(„) 
needs to be addressed due to its potential application to other operators such as those which 
contain four quark fields and one or more 7^-matrices. One reason for this is that the algebra of 
the r(„)-matrices has been widely studied, |1C, 11, 12, 13|, and also provides a simpler way to 
programme in a symbolic manipulation language than say the split 7^ algebra of [0]. Indeed 
the renormalization of j(„) at three loops which we carry out here is performed with the use of 
the Mincer package, [15|, written in the language Form, [16|. While the above motivations for 
this work have been in relation to QCD, we note that knowledge of the anomalous dimension 
of the currents j(„) are necessary for other problems. For instance, the relation between QCD 



and the non-abelian Thirring model, (NATM), in d-dimensions has been developed in [17| based 
on the earlier observations of |l^. To understand the connection further, one needs not only 
to have knowledge of the fundamental renormalization group functions such as the /^-function, 
but also the anomalous dimensions of the operators in both theories as well. By providing these 
for the n-dependent operators additional information is being determined for this area. For 
example, in the strictly two dimensional NATM the currents with n = 3 and 4 would be zero. 



However, in the d-dimensional context of the fixed point equivalence with QCD, [17, |^], they 
would be evanescent and important for the connection since from the QCD point of view they 
do not vanish in the limit to four dimensions. 

We now turn to the discussion of the computation. The procedure to renormalize the currents 
•?(n) ^" is to insert the operator in the quark two point function, use dimensional regularization 
and determine the poles in e, where d = 4 — 2e, before minimally subtracting them. In particular 
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the set of Feynman diagrams to three loops are generated with the Qgraf package, [19|, and 
converted into Form input for processing with the Mincer procedures, |15|. As this is a 
weU documented method we merely note the major issues in relation to the operators we are 
interested in. First, by constructing currents with the generalized 7-matrices F^^j"'^" defined by 



^(n) 



■ r 



(1) 



there is no mixing under renormalization between any of the Next to apply the MiNCER 

routines the Green's function must be converted to scalar integrals. This requires multiplying the 
quark two point function by F^^^^^ taking the spinor trace and dividing by the normalization 
tr(F^^j"'^"F(„)^j ^^). As we are only interested in flavour non-singlet currents for the moment, 
the 7-matrix strings contain either two F^^^'s or none. For those with none the trace operations 
is readily performed. For the former strings one has to use the properties of the F^^j"'^" which 
are given in [|l|, It transpires that we only need a subset of those results. By considering 

the Feynman diagrams which arise at three loops it is clear that the two F(„)'s are separated by 
at most six ordinary 7-matrices. Therefore, from the general expression, [12|, 



-p/il . . ./in . .Mm p 

(n) (m) (n)/ii.../i„ 



/(n,m)F-y- 



it is easy to deduce 

rf^)-'^"7^^r(„)^,...^„ = /(n,2)7^^ + (/(n,0)-/(n,2))r?^^ 

and 



(3) 



■^^n 



/(n,4)7'^7^V 
+ (/(n, 2) - /(n, 4)) [r/^''7'^7^ - 77'"^7'^7^ + ,^^7^7'^ 

+ (/(n, 4) - 2/(n, 2) + /(n, 0)) [r/^^r/'^^ - ^'^r/'^'' + 77^^r/'^'^] . 

(4) 

The analogous relation involving six 7-matrices between two contracted F(„)'s is too large to 
quote here though it clearly will involve the function /(n, 6). The values of f{n,m) we require 
are deduced from the general formula given in |12|. We find 
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They have been written in terms of /(n, 0) since this also occurs in the normahzing factor and 
hence wih cancel in the computation. With these lemmas the Tj-^-j's can be removed from all 
7-strings and the spinor trace evaluated in the normal fashion which then means the integrals 
are in the correct format for applying the Mincer procedures. 



Taking into account the three loop wave function renormalization, |2C], of the quark fields 
present in the original operator, we find the following result for the anomalous dimension of 
for arbitrary n, 

n — l){n — 3)CFa 

4(n - 15)TFNf + (ISn^ - 12671^ + 163n + 291)Ca 

(n - l)C7^a2 



7(n)( 



+ 
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15846n^ + 15933n + 11413 



- 9(n - 3)(5n^ - 20n + 1)Cf 
(l44n^ - 158471"^ + 6810ri 

- 216n(n - 3)(n - 4)(2n2 - 8n + 13)C(3)) C 
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432n(n-3)(n-4)(3n2 
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- 12n + 19)C(3)) CaCf 
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+ O(a^) 
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where a = as/47r = 51^/(1677^) is the coupling constant, Tp, Cf and Ca are the usual colour 
group Casimirs and is the number of quark flavours. There are various checks on this result. 
First, as the operator itself is gauge invariant, its anomalous dimension must be independent of 
the covariant gauge parameter, ^. Therefore, in the calculation we used a gluon propagator of 
the form (r/^,y — ^p^Pu/p^) /p^ where p is the momentum, and observed that ^ cancelled in the 
final result, (^). Second, the known three loop MS results for n = 0, |§ ^, and 1 emerge. In 
the latter case the anomalous dimension vanishes as the current corresponds to the conserved 
electric current. Moreover, the value of (^) for n = 3 and 4 corresponds to the naive values 



obtained by Larin in 
current ^a^^ip is 

7(2) (a) = 



Hence, we can deduce that the anomalous dimension for the tensor 



Cpa + [257Ca 
13639C1 - 



niCp - 52TFNf] 



18 



+ 



4320C(3)Cl + 12096C(3)CaCf 
20469CaCf - 1728C(3)CArpiV/ - AQlQCATpNf 
6912C(3)C|^ + 6570C| + 1728C(3)CFTpiV/ 

Cpa^ 



108 



+ 0(a^ 



+ mQCFTFNf 

by substituting n = 2 in (^). The first two terms are in agreement with 
group is SU (3) we have 

4a 2[26iV/ - 543] 

y 27 



(7) 

When the gauge 



7(2) (a) 
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367V/ + 1440C(3)iV/ + 5240iV/ + 2784C(3) - 52555 — + O(a^) (8) 



a3 



or numerically 



72(a) = 1.333333a - (1.925926iV/ - 40.222222)0^ 

_ (0.444444A/ + 86.061259iV/ - 607.512020)0^ + O(a^) . (9) 

We note that whilst we have concentrated on the flavour non-singlet current, the expression 
d?]) also corresponds to the anomalous dimension of the (anomaly free) flavour singlet tensor 
current. This follows trivially since the Feynman diagrams with a closed quark loop with one 
ipa^'^'ip operator insertion in it, has an odd number of 7-matrices and are therefore zero upon 
taking the spinor trace. 

One question which naturally arises out of choosing to compute these anomalous dimensions 
with the generalized 7-matrices, r(„), is that of how they relate to currents which involve 7^. As 
is well known its treatment in dimensional regularization is a technically difficult exercise due to 
the fact that it has no natural d-dimensional generalization, [|l^]. In this issue of renormalizing 
the axial vector and pseudoscalar currents within the automatic multiloop computer algebra 
approach was addressed. To treat such currents one first of all defines 7^ in terms of e^iye^pT'^^^'' 
where e^^ap is the totally antisymmetric four dimensional pseudotensor, Then products of the 
e-tensor are replaced by a function of the metric, rj^i,, which does have a natural d-dimensional 
extension. To compensate for the lack of continuity of the 7^ definition in ci-dimensions a finite 
renormalization is performed in addition to the usual MS subtraction of the Green's function, 
1^]. The condition for the finite renormalization is to impose the obvious anticommutativity 
of 7^ on the finite renormalized Green's functions. We have summarized the renormalization 
of these quark currents of previous work, to allow the interested reader to contrast their 
renormalization here. Clearly within the context of the complete set of currents j^^^"^" the 
correct anomalous dimension for the axial vector and pseudoscalar currents should somehow be 
present in our results. We now address this. 

In choosing to work in a spacetime which is d-dimensional, with d non-integer, the problem 
of defining 7^ is in some sense bypassed in that one can regard it as an object which does 
not exist naturally. Also the ordinary 7-matrix basis ceases to be finite dimensional. Indeed 
in this situation the r(„)-matrices provide a more natural basis for performing d-dimensional 
calculations. (For a more detailed discussion on this point see, for example, [|^, 13, ^.) As 



7 is a manifestly four dimensional object one need only be concerned about incorporating its 
effect when projecting from the d-dimensional spacetime, with its infinite F^^^-basis, onto the 
integer dimensional physical spacetime. For us this projection is not achieved by analytically 
continuing the product e^^ape"^^^^ to d-dimensional spacetime, [^J. Instead within the context 
of the infinite dimensional 7-matrix basis one determines a finite renormalization constant Z^^^ 
from a condition similar to that of [^], 

G(2)(p2,n) = Zfi,G(2)(/,4-n) (10) 



where G^^^ {p^ ,n) is the finite part of the quark two point Green's function after minimal subtrac- 
tion where the current has been inserted. In other words it is the renormalized Green's 
function. Also in ( p!o[) it is understood that all the contributing diagrams have been multiplied 
by [n) ^i]_...^in before taking the spinor trace and dividing by the appropriate normalization men- 
tioned earlier. Choosing the argument of the right side of (|lO|) to be (4 — n) ensures that for 
instance the finite renormalization constant which emerges for the n = 4 current will give the 
same anomalous dimension as the n = current with no finite renormalization. Therefore from 
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our calculations we find that 



+ 



1 + A{n-2)CFa 

(36n2 - 144n + 1)Ca 

216(6n^ 

- 3(144n^ 
+ 216(6n^ 

+ 6(108n^ - 1080n^ + 4169n^ - 6314n + 200)CaCf 

- A32{9n^ - 72n^ + 200n2 - 224n + 57)C(3)CaCf 

- 8{6n^ -24n + W7)CATFNf - l728C{3)CATFNf 
+ 54(17n^ - 76n^ - 32n^ + 352n - 76)C' 



18n(5n - 16)Cf + 4TFiV/ 

48n3 + 134n2 - 152n + 39)C(3)Cl 

1152n^ + 3292n2 - 3952n - 479)Cl 
48n^ + 134^2 - 152n + 39)C(3)C7^ 
1080n^ + 4169n2 
72n^ + 2mr? 



(n - 2)CFa^ 
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■^2 



+ 2592(n^ 
+ 12(34n2 



8r? + 22n2 - 24n + 6)C(3)C| 



148n + 145)Cp^rpiV/ 



+ 1728C(3)CFrFA^/ 



208T|.7V/ 



(n - 2)Cfc? 
81 



(11) 



and we note that it is independent of the gauge fixing parameter. Moreover, by construction 
it evaluates to unity for n = 2 as it ought. Hence, to deduce the anomalous dimension of the 
currents '07'^7^^ and '^'~f'^ themselves by this method, one must therefore add the piece 



dln^fin 



4(7f(11Ca - 4rF iV/)(n - 2)a? 



+ 2Cf 



(396^2 



3 

- 1584n 



198(5n2 - 20n + 8)CfCa 



+ 72(5n2 



- 16(9n^ - 36n - 
20n + ll)CFTFNf 



2b)CATFNf 
{n - 2)a^ 
27 



(12) 



to 7(3) (a) and 7(4) (a) respectively for n = 3 and 4 where /i is the renormalization scale. This 
effectively projects out the true component of (^) for four dimensional spacetime and correctly 
reproduces the known results of 0, |5|. We note that the three loop contribution of (12) only 
involves the two loop term of (pT|). Also if one was working with reference to two dimensions 
the argument of the right side of (|l^) would instead be (2 



n 



Finally, we briefy comment on this approach to flavour singlet current anomalous dimen- 
sions which we have also analyzed. As noted earlier the flavour singlet and non-singlet currents 
coincide for even n. However, if one were to use the r(„)-basis to study the anomalous di- 
mension of singlet currents for n odd then there are two further issues to be dealt with. The 
first is the evaluation of 7-strings with only one r(„)-matrix in the spinor trace. To handle 
these diagrams the following results are necessary. First, one decomposes the 7-string arising 
from the propagator and vertices into the F(„)-basis and then uses the general property that 



tr(F 



(n) 



(P) 



ti'(F(„)^^...^„F' 



(Tl. 

(9) 



is proportional to Spq. The tensor of proportionality in- 



volves a linear combination of products of {p + q)/2 ry-tensors which respect the antisymmetry of 
the Lorentz indices of F(p) and F(q). At three loops the largest value of p or (7 is 5 which means 
that a small set of traces needs to be deduced. This is achieved by noting that, 13 1, 



tr frf 



Y.Ul..Mr 



^ = 4(-l)"(" ^^^'^6mnnl['n^^''^ . ■ ■rj^""" + antisymmetric permutations] 



(13) 
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With these identities the graphs with a single operator insertion can be determined and the 
basic anomalous dimension analogous to (^) deduced as 



sinelet / \ 



7(„)(a) + l2dn,3CFTFNfa 
218 

' 3 



+ 



Sn,3CF ( ^CaTfNj + ^TliV/ - 36Ci.ri.iVj 



+ <5n,5(480C(3) + 80)- 



^fund 



where 



df^ 



Tr ( T^°-T^T'^^ 



(14) 



(15) 



are the colour group generators and iVfyjj^j is the dimension of its fundamental representa- 
tion. The 5-symbols arise from the spinor traces with one r(„)-matrix. Whilst an additional 
colour group Casimir enters for the case n = 5, this current is evanescent with respect to four 
dimensions and is therefore not physically important. Clearly to relate the n = 3 current to 
the four dimensional axial vector current anomalous dimension one requires an additional finite 
renormalization and this arises from two places. The first occurs by imposing the condition ([lO| ) 
and deserves little more comment only to record that an extra term arises in (^) when n = 3 
which is 



2Ci. 



IICaTfNj 



28r|iV/ 



(16) 



However, this additional finite renormalization will not ensure the correct singlet anomalous 
dimension emerges since the (four dimensional) chiral anomaly will not be preserved. This 
requires another finite renormalization and is achieved here by the method developed in Q but 
with the anomaly equation treated inside a quark two point function in contrast to the gluon two 
point function considered in [^. As this procedure only involves the gluonic operator inserted 
in a one loop diagram to the order we are interested in, we note that the additional contribution 
will be 



2Ci. 



+ 



2AT],Nf 



(17) 



66CATi.iV/ 

which gives the correct singlet axial vector current anomalous dimension, ^. 

In conclusion we have provided a new term in the series for the anomalous dimension of 
the tensor current in QCD in the MS scheme. Also by considering the generalized 7-matrix 
basis, r^^), we have demonstrated how the anomalous dimensions of the currents which involve 
7^ emerge. Indeed we believe this is an important aspect of the calculation as it in principle 
provides a more systematic and alternative way of renormalizing currents or other composite 
operators in QCD which involve the purely four dimensional object 7^. For example, one need 
only calculate Green's functions with a general insertion and then the anomalous dimensions 
for a variety of operators will emerge by choosing the appropriate value of the parameter of 
the inserted operator. It would be interesting, for instance, to develop this approach for more 
complicated and physically important operators such as the four quark operators which are 
fundamental to evaluating QCD corrections to weak processes. 
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